Inflation-Produced Magnetic Fields in Nonlinear Electrodynamics 
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We study the generation of primeval magnetic fields during inflation era in nonlinear theories of 
electrodynamics. Although the intensity of the produced fields strongly depends on characteristics 
of inflation and on the form of electromagnetic Lagrangian, our results do not exclude the possibility 
that these fields could be astrophysically interesting. 
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I. INTRODUCTION 



All galaxies seem to be permeated by magnetic fields 
with intensities of order -B ga iactic — 10 _6 G 

To explain the galactic magnetism, generally one needs 
the presence of seed magnetic fields prior to protogalaxy 
collapse. When a protogalaxy collapses to form a galac- 
tic disk, magnetic fields suffer an amplification (mainly 
due to magnetic flux conservation) of order „4 pg ~ 10 4 
Q. Moreover, due to magnetohydrodynamic turbu- 
lence effects and differential rotation of galaxy, seed 
fields can be further amplified. This last mechanism, 
know as "galactic dynamo" 0], can be very efficient 
and, in principle, allow extremely week seeds to repro- 
duce the properties of presently-observed galactic fields. 
Dynamo action produces an exponential amplification 
Adyn — e rAt , where the growth rate T is a model- 
dependent quantity, and At = tf — U is the time dur- 
ing which the dynamo operates. The minimum and 
maximum values of Y that can be found in the litera- 
ture are V ~ O^Gyr" 1 (H and V ~ SGyr" 1 |. In a 
spatially-flat Friedman universe, the time interval is Q 

At 
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Here, Hq = 100/ikmsec 4 Mpc is the Hubble constant 
(h ~ 0.73 0), tj u = (1 + z 4j )- 3 / 2 (fVan) 1/2 , z is 
the red-shift, and f2 m ~ 0.28 and Q\ ~ 0.72 [7[ are the 
actual energy densities, in units of the critical density 
p CI = 3Hq/8ttG, associated to matter and cosmological 
constant. Un upper bound on the value of red-shift at 
which the dynamo begins to operate, Zi, is given by the 
red-shift z pg ~ 50 at which a protogalaxy separates from 
the Hubble flow to then collapse [6]. If dynamo is effi- 
cient during galactic disk formation or if its amplification 
becomes effective only after that, is still not clear. In the 
latter case, since astronomical observations indicate that 
disk galaxies at z ~ 3 are still in progress of being formed 
3, one should take this value of red-shift as a conserva- 
tive bound on Z{. A lower bound on zj is z/ ~ 0.4, since 
microgauss magnetic fields have been detected in galaxies 
at that red-shift [§]. 



The galactic magnetism can be then explained as the 
result of the amplification of comoving seed fields as 



> 10" b (l 



z ps) -^pg Aiyn ' where 
pg )~ 2 takes into account the adiabatic 



strong as B f 
the factor (1 + z 
scaling of the magnetic field from the protogalaxy col- 
lapse until today. Taking (T, Zi, Zf) = (5, 50, 0.4) we have 
B sccd > 10~ 33 G, while for {T,Zi,z f ) = (0.45,3,0.4) we 
get -Bsced ^ 10 _15 G. In order to have an efficient galactic 
dynamo, however, the seed magnetic field must be corre- 
lated on comoving scales of order lOkpc. 
We observe, also, that without dynamo amplification 
a comoving seed field as strong as £? scc d ^ 10~ 14 G is 
needed to explain galactic magnetism. In this case, how- 
ever, the field must be correlated on comoving scales of 
order of linear dimensions of a protogalaxy, that is IMpc. 

Surprisingly, nanogauss magnetic fields correlated on 
scales of order IMpc, have been also detected in galaxy 
clusters and superclusters. This last observation seems to 
indicate that the entire universe is magnetized 0, [Tol . 11 1 . 



Essentially, there are two possible classes of mecha- 
nisms to produce cosmic fields depending on when they 
are generated [l2l | : Astrophysical mechanisms acting dur- 
ing large-scale structure formation 13j , and mechanisms 
acting in the primordial universe, during 0E1EEE1 
18 . 13 . 20l . 2lf or before [22] inflation. However, we can 



admit the existence of strong fields in the primordial uni- 
verse provided that their presence do not spoil predic- 
tions of the standard cosmological model, such as that 
of Big Bang Nucleosynthesis (BBN) [23j, Large-Scale- 
Structure Formation (LSS) j^, and Cosmic Microwave 
Background (CMB) [25| . It turns out that limits coming 
from LSS and CMB are more stringent than those from 
BBN. Putting together the limits found in Refs. (23,125 1 
it results that, for comoving scales in the range 400pc < 
A < d H (t ), where d H {t ) ~ H^ 1 - 4000Mpc is the 
present size of the universe, the maximum strength al- 
lowed to comoving primordial fields is B ~ 10 _9 G. 

In the ambit of generation of cosmological fields in the 
early universe, the mechanisms operating during infla- 
tion are particularly attractive since they produce large- 
scale correlated fields. Magnetic fields created after in- 
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nation, instead, suffer from a "small-scale problem" , that 
is their comoving correlation length is much smaller then 
the characteristic scale of the observed cosmic fields [how- 
ever, if magnetohydrodynamic turbulence operates dur- 
ing their evolution, an enhancement of correlation length 
can occur (see, e.g., Ref. |26l|)]. 

It is worth noting that, due to conformal invariance of 
standard (Maxwell) electrodynamics and to the fact the 
spacetime described by the Robertson- Walker metric is 
conformally flat, magnetic fields generated at inflation 
are vanishingly small. For this reason, all generating 
models proposed in the literature repose on the break- 
ing of conformal invariance of Maxwell theory. This has 
been attained, for instance, by non-minimally coupling 
the photon with gravity [l_4], introducing interactions of 
photons with scalar, pseudoscalar, or vector fields (such 
as inflaton Qjl, dilaton fltjj| . pseudo-Goldstone bosons 
17 1, axion or "graviphoton" [lj|), taking into ac- 
count the so-called Quantum Conformal Anomaly [20l l . 
and so on 



21] 



In this paper, we study the possibility to generate seed 
magnetic fields during inflation in nonlinear theories of 
electrodynamics (NLE) described by the general action 



S 



1 

47T 



(1) 



where F = jF^F^ ', with F, 



d^A^ the elec- 



fi A 

tromagnetic field strength tensor, and g — det ||<7p„|| is 
the determinant of the metric tensor. 
Since the standard Maxwell theory is a good theory for 
low energies, we shall assume that nonlinear Lagrangians 
reduce to the Maxwell one, C(F) ~ —F, in the limit of 
small fields F. 

A considerable amount of interest has emerged in the 
last few years in cosmological effects of nonlinear elec- 



trodynamics [27J, [28| . This is due principally to the fact 
that some theories of NLE are able to produce inflation 
[H, [2§| , a period of cosmic acceleration [28|, [3(1 , and can 
also avoid the problem of initial singularity |28l . 31 , 32| • 
In general, nonlinear electrodynamic theories are non- 
conformally invariant. As we shall see in the next Sec- 
tion, depending on the actual form of the Lagrangian, 
astrophysically interesting magnetic fields can be gener- 
ated during inflation. 



II. GENERATION OF SEED FIELDS IN NLE 
A. Equations of Motion 

We will work in a flat universe described by a 
Robertson- Walker metric, ds 2 ~ a 2 (r])(drj 2 — dx 2 ), where 
a(j]) is the expansion parameter and rj is the conformal 
time related to the cosmic time t through dr/ = dt/a. 
Introducing the electric and magnetic fields E and B in 



indices range from 1 to 3), and varying the action (JTJ) 
with respect to A^, we get the equations of motion: 



<9(a 2 E) 
di] 



- V x (a 2 B) = - 



2T1 ,_ 9 l n ^£ a 2 E + Vln£FXa 2 B) 
OX] 

(2) 

and V • (£j?E) = 0, together with the Bianchi identities: 
d v (a 2 B) + V x (a 2 E) = and V • B = 0. Here, subscript 
on C denotes differentiation, and spatial derivatives are 
taken with respect to comoving coordinates. 
We are interested to the evolution of electromagnetic 
fields outside the horizon, that is to modes whose phys- 
ical wavelength is much greater then the Hubble ra- 
dius <iff, Aphys > dn, where A p h ys = a\, dn ~ H~ l , 
and A is the comoving wavelength. Since rja ~ H , 
introducing the comoving wavenumber k = 2tt/X, the 
above condition reads \krj\ <C 1. Observing that the 
first Bianchi identity gives B(fcry) ~ krj^kr]), we have 
that B 2 is negligible with respect to E 2 , and we can 
write d v \n£ F ~ — (<9, ; E 2 /2) (d\n£ F /dF). Moreover, we 
can neglect the second term with respect to the first 
one both in the left- and right-hand-side of Eq. ([2]). In 
fact, it results |V x a 2 B|/|5 I? (a 2 E)| ~ \kr)\ 2 < 1, and 
|Vln£ F x a 2 B|/|(cyn£ F )a 2 E| ~ \kr]\ 2 < 1. Finally, 
multiplying Eq. ([2]) by E, and solving with respect to 
F ~ -5E 2 , we get 



(C F ) 2 F cx a 



(3) 



Knowing the form of nonlinear Lagrangian (see below), 
from the above equation one gets the evolution law for 
the electric field outside the horizon. Consequently, us- 
ing the first Bianchi identity, one finds how super-horizon 
magnetic fields scale in time during inflation (see subsec- 
tion D). 



B. Initial Electromagnetic Spectrum 

During inflation, all fields are quantum mechanically 
exited. Because the wavelength A associated to a given 
fluctuation grows faster then the horizon, there will be a 
time, say t\, when this mode crosses outside the horizon 
itself. After that, this fluctuation cannot collapse back 
into the vacuum being not causally self-correlated, and 
then 



the usual way as Foi 



-a 2 Ei 



e ljk a 2 B k (Latin 



''survives" as a classical real object [33j. The en- 
ergy associated to a given fluctuation is subjected to the 
uncertainty relation, AEAt > 1. Therefore, the energy 
density in the volume AV, £ — AE/AV, is approxima- 
tive^ given by £ ~ H A , where H is the Hubble parame- 
ter. Here, we used the fact that at the horizon crossing 
At ~ H- 1 and AV - H~ 3 When a comoving 

length A crosses the horizon it results \kr/\ ~ 1, and then 
from the first Bianchi identity we get B 2 (A) ~ E 2 (A). 
Therefore, since F = — ^(E 2 — B 2 ), at the crossing the 
nonlinear terms in the electromagnetic Lagrangian are 
negligible, and the energy density is simply given by 



3 



E ~ (1/8tt) (E 2 



Consequently, the spectra of the 



quantum-mechanically generated electric and magnetic 
fluctuations are given by 



B 2 (A) 
4tt 1 



E 2 (A) 
4tt 



Pt 2 ot(A) 



liii 



(4) 



'pi 



where in the last equality we used the Friedman equa- 
tion H 2 = (8ir / 3) pt t / rrip V Here, ptot is the total energy 
density during inflation and mpi ~ 10 19 GeV the Planck 
mass. 

In the following we shall consider both the case of de Sit- 
ter inflation and the case of "Power-Law Inflation". It 
is useful to write the expansion parameter as a(i]) oc if . 
During de Sitter inflation s = — 1 and the total energy is a 
constant, 1 ptot = M 4 . For power-law inflation described 
by the equation of state ptot = 7Ptot with —1 < 7 < 1/3, 
we have s = 2/(1 + 37) and the total energy decreases 
as a power of the expansion parameter, ptot oc a -3 ' 1+7 -*. 
Therefore, the total energy, when a comoving length A 
crosses outside the horizon, depends on A [see Eq. (TIB"]) ]. 



C. Plasma Effects 

The effects of a conducting plasma in the early uni- 
verse on the evolution of magnetic fields, are taken into 
account by adding to the electromagnetic Lagrangian the 
source term JM M [14| . Here, the external current J M , 
expressed in terms of the electric field, has the form 
J M = (0, (T C E), where a c is the conductivity. Plasma 
effects introduce, in the left-hand-side of Eq. (|2|), the 
extra term — acx c (a 2 E). In the limit of high conductiv- 
ity, a c — > 00, one finds E — ► and, consequently, from 
the first Bianchi identity, it follows that the magnetic 
field is frozen into the plasma and evolves adiabatically, 
a 2 B ~ const. More precisely, the solution of Eq. ([2]), for 
modes outside the horizon and when plasma effects are 



taken into account, is (Cf) F 



oc a 



1 exp(— 2 J di] aa c 



Observing that Jdr]aa c ~ a c /H, we get that the mag- 
netic field can be considered as frozen into the plasma 
when a c H . After inflation, the universe enters in the 
so-called reheating phase, during which the energy of the 
inflaton is converted into ordinary matter. In this paper, 
we shall restrict our analysis to the case of instantaneous 
reheating, that is after inflation the universe enters the 
radiation dominated era. In this era, the conductivity is 
approximatively equal to a c ~ T [14j . while H ~ T 2 /mp\ 
[35|] , where T is the temperature. Hence, the condition 
of freezing of the magnetic field becomes T <C mpi. 



The spectrum of gravitational waves generated at infla- 
tion is submitted to constraints coming from CMB analy- 
sis which requires ptot(A) to be less than about 10 _8 m pl 
on the scale of the present Hubble radius [3]. This, 
in turns, converts in a upper limit on the value of M, 
M < 10~ 2 topi. (One must impose also that M> IGeV, 
so that the predictions of BBN are not spoiled [bj].) Since 
the temperature at the end of inflation is T en( j = M, we 
conclude that after inflation the universe is a good con- 
ductor, er c 3> H, and the magnetic field evolves adiabat- 
ically, irrespective of when it (eventually) reenters the 
horizon. 



D. Form of NLE Lagrangian 

In this paragraph, we consider three models of non- 
linear electromagnetic Lagrangian. In all cases, the La- 
grangian depends on a free mass parameter, m, such that 
in the formal limit m — > 00 we recover the standard 
Maxwell theory. More precisely, it results C(F) ~ — F 
for \F\ -C ra A . In the case of small fields, |.F| < m 4 , 
inflation-produced fields are vanishingly small, and then 
cannot explain the presently-observed fields. For this rea- 
son, we shall restrict our analysis to the case of strong 
fields, |F| > m. 

As a first model, we consider the family of Lagrangians 



C(F) 



-F 



(5) 



i=2 



where i takes values on the integers, and the coefficients 
Cj have dimension [Mass] 4 ^ 1- -'-'. We assume that Cj = 
m^^'dj, where dj are dimensionless constants of order 
unity. In a cosmological context, this type of Lagrangian 
for n = 2 has been widely studied in the literature (see, 
e.g., Ref. [l^,!!^]). In Ref. [36j |. it has been shown that 
the Lagrangian 2 



C KK (F) = -F + T (b-l)F 2 



(6) 



T being a parameter with dimension [Mass] -4 and b a 
dimensionless parameter, derives from higher-curvature 
gravity in Kaluza-Klein theory. 

In the limit of strong electromagnetic fields, we have 
C(F) ~ c n F n . In this case, Eq. © gives 



E 



-4/(2n-l) 



(7) 



1 In de Sitter inflation, the spectrum of electromagnetic fluctua- 
tions when crossing the horizon is ~ \F I1V \/H ~ H, that is a 
scale-invariant spectrum corresponding to the Gibbons-Hawking 
temperature Tqjj = H/(2n) |34| . 



2 The full Lagrangian is C KK (F, G) = -F + T[(b-1)F 2 - 3G 2 /2], 
where G = \F^F^ V = E-B, and F^ = (1/2^=9 ) <lV- v f>° F pa is 
the dual electromagnetic field strength tensor. However, in this 
paper, we are concerned only with nonlinear theories depending 
on the invariant F. 



4 



where a\ — a(ii), and as initial value for the electric field 
we have taken that at the horizon crossing, Ei = E 2 | tl . 
From the first Bianchi identity we get a 2 B ~ (krj)a 2 E + 
const, while from Eq. ([7]) we have a 2 E ~ (krj)£ with 
£ = 4s(n— l)/(2n— 1). Since we are assuming \kr)\ -C 1, 
if £ < — 1 the magnetic field evolves as 

/ >. 2(2n-l-2s)/s(2n-l) 

b^b;(£) ( s) 

while, if £ > —1, it scales adiabatically. For n — 1 we 
have £ = while for n > 2 it results £ < — 1. Moreover, 
since 2(2n — 1 — 2s)/s(2n — 1) > —4 for n > 2, a "supera- 
diabatic amplification" {i.e. B 2 evolves less slowly than 
the usual a -4 ) occurs during inflation. 

We now consider a "toy model" described by the La- 
grangian 



C(F) = -Fe 



-cF 



(9) 



where c : 



l d, with d > a dimensionless constant of 



order unity. The exponential self-coupling in Lagrangian 
© resembles to the exponential coupling oc F^F^e "^ , 
a being a dimensional constant, between the inflaton 
(dilaton) <f> and the electromagnetic field, introduced in 
Ref. [15| (Ref. |l6|). In our case, the scalar field is re- 
placed by the scalar F. 

In the limit of strong fields, the solution of Eq. © is 
approximatively given by E 2 ~ E 2 + to 4 \n(ai/a) 4 ^ d . Ne- 
glecting the logarithmic term, we have that the second 
model is equivalent to the first one with n — > oo. 

As a third model, we consider the Born-Infcld (BI) 
Lagrangian 3 



£ BI (F) =m 4 1 




(10) 



where, for all field configurations, the condition 2F/m 4 > 
— 1 has to be satisfied. Born and Infeld proposed their 
theory j3?| in order to eliminate the divergence in the 
energy of a point-charge particle. Indeed, in this the- 
ory the self-energy of a point-like charge is always finite 
and proportional to the parameter m. Born-Infcld model 
also appears in quantized string theory (38j (in this case, 
rn 2 = 27ra', where a' is the string tension parameter). 
Cosmological effects of BI electrodynamics, such as gen- 
eration of an inflationary phase [29], have been deeply 
studied in recent years. 

Equation © gives E 2 = m 4 [(a/ai) 4 (m 4 /E 2 - 1) + l]" 1 . 

then for a > ai we get the usual behav- 



If E 2 < to 4 
ior E 2 cx a~ A . 



If Ef 



we have E 



i 4 for all 



times. Therefore, this model is equivalent to the first 



one with 



oo, and with the condition \F\ > 



that 



is Ej > to 4 , replaced by E 1 



3 The full Lagrangian is C m (F,G) = m 4 [1 - ^l + 2F/m 4 -G 2 /m s ] 
(see also footnote 2). 



E. Present Magnetic Fields 

We now derive the actual strength of magnetic fields 
generated during inflation in nonlinear theories described 
by the above three Lagrangians. 

During de Sitter inflation any weak field, \F\ < to , 
is exponentially washed out. We then analyze the case 
in which electric fields remain strong from the first hori- 
zon crossing to the end of inflation. Observing that the 
electric field is a non-increasing function of time, we then 
assume that the quantity E 2 (A)/TO 4 | 4end is greater then 1, 
where i e nd is the time corresponding to the end of infla- 
tion. Taking into account Eqs. (J4| and ([7]), we can write 
the above condition as 



/ M \ 2 ~ M 

to < 10 20 (10 24 A 10kpc )^ GeV, (11) 

V TOPI 



where \i = l/(2n — 1), and we used a cn d/ai = e N ^ ~ 
10 24 Aiokpc-^/™pi, N(\) being the number of e-folds 
elapsing from the crossing of a comoving length A outside 
the horizon to the end of inflation 



1$. 



The actual value of the magnetic field follows from Eq. ([8|) 
and is -Bnle — -Bie _/3A,( ' A '(a cn d/ao) 2 where B±, the mag- 
netic field strength at the horizon crossing, is given by 
Eq. ©, and j3 — 1 + 2/i. The last term in the above equa- 
tion takes into account the adiabatic dilution of the mag- 
netic field from the end of inflation until today, a = olq. 
Using the relation (valid during radiation and matter 
dominated eras) a oc g^' 3 T _1 , g*s(T) being the num- 
ber of effectively massless degrees of freedom referring to 
the entropy density of the universe (3o| , we arrive to 



B 



NLE 



10- 4 (lO^A^kpc)-' 3 



(-) 



2-/3 



G, (12) 



where we used the values To ~ 2.35 x 10 13 GeV, 
<7*s(To) ~ 3.91, and we assumed #*s(T c „ d ) - 10 2 j^. (It 
is useful to know that 1G ~ 6.9 x 10~ 20 GeV 2 .) Observe 
that the magnetic field during de Sitter inflation evolves 
as B oc a^' 3 . The case of standard electromagnetic La- 
grangian corresponds to (3 = 2. Therefore, one finds 
for A = lOkpc the vanishingly small value B ~ 10 _52 G. 
In nonlinear electrodynamics, taking M ~ 10~ 2 topi and 
A = lOkpc, we find that -Bnle is an increasing func- 
tion of n. In particular, we get Bnle — 10 _45 G for 
n = 2, B NLE ~ 10- 33 G for n = 8, and B NLE ~ 10- 30 G 
for n — > oo, together with the conditions to < 10 8 GeV, 
to < 10 14 GeV, and to < 10 16 GeV, respectively. 
From the above results, we see that NLE effects are able, 
in principle, to produce magnetic fields that can seed 
galactic dynamo. 

In the case of power-law inflation, a comoving length 
A crosses outside the horizon when [14| 



p tot (A)| tl ~ ( 10 24 — A 10kp c 

TOPI 



-2.r 



M 4 , 



(13) 
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TABLE I: Actual strength of the magnetic field, Bnle, pro- 
duced during power-law inflation via nonlinear electrody- 
namic effects at the comoving scale A = lOkpc. The last 
two rows refer to A = IMpc. The index n and the mass m de- 
fine a particular choice of the nonlinear Lagrangian (see text 
for details), while M 4 is the total energy density at the end of 
inflation. The magnetic field depends on the comoving scale 
A as -Bnle(A) oc A _/3 . 



n 


M(GeV) 


Bnle(G) 


m(GeV) < 





2 


10 16 


1Q -43 


10 s 


1.6 


2 


10 12 


lO" 37 


10 s 


1.1 


2 


10 9 


1Q -32 


10 3 


0.6 


3 


10 16 


lO" 37 


10 11 


1.3 


3 


10 12 


1Q -30 


10 s 


0.7 


3 


10 9 


1Q -26 


10 6 


0.1 


3 


10 6 


10~ 23 


10 3 


-0.7 


00 


10 16 


1Q -28 


10 15 


0.9 


00 


10 12 


1Q -20 


10 13 


0.1 


00 


10 9 


1Q -14 


10 13 


-0.6 


00 


10 8 


10" 12 


10 12 


-0.9 



where M 4 = ptot(A)| en d is the total energy density at the 
end of inflation, and x = 3(1 + 7)/(l + 37). The bound 
on graviton production previously discussed translates to 
[3] x > a; min , where x min ~ [4+21og 10 (M/m P i)]/[29.8 + 
log 10 (M/mpi)]. In the following, we assume for simplic- 
ity that x = iCmin, and that the electric field remains 
strong during inflation. This corresponds to the "best 
case scenario" , or to the minimum dilution of the mag- 
netic field during inflation. Taking into account Eqs. ((3]) 
and 0, and using a cnd / ai = (M* / p tot (X)\ tl )~ 3( - 1+ ^ , the 
condition E 2 (A) /m 4 1 tond > 1 translates into Eq. (jlip with 
the replacement /u —* // = fi + (1 — /i)a; m m- The actual 
strength of the inflation-produced magnetic field follows 
from Eqs. ©, (JTHJ) , and is given by Eq. (JUJ with 
(3 replaced by (3' = 1 + 2//. The standard case of linear 
electrodynamics corresponds to /?' = 2, and then reduces 
to that studied for de Sitter inflation. 
In Table 1, we show the values of -Bnle (for the case 
of power-law inflation) for different values of n and M, 
together with the condition on m in the order that elec- 
tromagnetic fields be strong. Looking at the Table, we 
see that magnetic fields able to seed galactic dynamo or 
directly explain galactic magnetism can be produced. 

III. CONCLUSIONS 

Large-scale magnetic fields are ubiquitous in the 
present universe. Astrophysical observations have proved 
the existence of microgauss magnetic fields in all types 
of galaxies (spiral, elliptical, barred and irregular). Re- 
markably, nanogauss fields have been detected in galaxy 



clusters, and probably in superclusters, with correlation 
lengths of ~ IMpc. This sort of "cosmic magnetism" 
could have been arouse out of quantum electromag- 
netic fluctuations excited during an inflationary epoch 
of the universe. However, in standard, conformally- 
invariant theory of electrodynamics, inflation-produced 
fields are vanishingly small, and then cannot explain the 
presently-observed fields. Nevertheless, a lot (sometimes 
exotic) mechanisms able to break conformal invariance of 
Maxwell's electrodynamics, and then to produce astro- 
physically interesting fields, have been proposed in the 
literature. 

In this paper, we have investigated the possibility to 
generate magnetic fields during inflation era in nonlinear 
theories of electrodynamics, in which conformal invari- 
ance is naturally broken. We have found that, for a wide 
range of parameter space of inflationary models, mag- 
netic fields of cosmological interest can be created. In 
particular, we have shown that magnetic fields able to 
seed galactic dynamo or to explain directly the galactic 
magnetism could be a natural consequence of such theo- 
ries. However, since our results strongly depend on the 
actual form of the (unknown) nonlinear electromagnetic 
Lagrangian, they cannot give a definitive answer to the 
question "Why is our universe magnetized?" 

L.C. would like to thank Z. Berezhiani, D. Comelli, 
M. Gasperini, E. Lisi, and F. L. Villante for useful dis- 
cussions. 

Note Added. After completion of this paper, we 
noted a recent preprint by K.E. Kunze [arXiv.0110. 24-351, 
who exploits a very similar idea, but using a differ- 
ent (power-law) parametrization for the nonlinear la- 
grangian. We agree with Kunze 's results in the common 
subcase of integer- exponent monomial. We note, how- 
ever, that we can obtain cosmic magnetic fields of the 
observed size via power-law inflation (see Tab. I and re- 
lated comments), without necessarily invoking a galactic 
dynamo as in Kunze 's paper. 
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